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Abstract. - We present exact expressions for the quantum sloshing of Bose-Einstein condensates 
in a tilted two- well potential. Tunneling is suppressed by a small potential difference between wells, 
or tilt. However, tunneling resonances occur for critical values of the tilt when the barrier is high. 
At resonance, tunneling times on the order of 10-f 00 ms are possible. Furthermore, such tilted res- 
onances lead to a dynamical scheme for creating few-body NOON-like macroscopic superposition 
states which are protected by the many body wavefunction against potential fluctuations. 



Bose-Einstein condensates (BECs) in optical lattices are 
^— ^ an ideal medium for studying a vast range of quantum 
I many-body phenomena [T][2] , including macroscopic quan- 
V^tum tunneling. A two- well potential is a simple limiting 
case which nevertheless exhibits rich quantum behavior. 
Spatially separate BECs in a two- well potential have been 
O c reated in experiments [3] and tunneling times on the or- 
^Jjder of 50 ms have been observed [3]. The lifetime of a 
£h typical experiment is 1-100 s [5]- This system, in certain 
^limits (see Eq. (JTJ) below), maps onto the Lipkin-Meshkov- 



Q-Glick (LMG) model [S[7]|ll the A-body generalization of 
• »the two-state problem, in which N particles can occupy 
^fwo single-particle modes. Dynamical instability at the 
level of mean field theory (MFT) is associated with gen- 
5-h e'ration of strongly correlated states on a microscopic ex- 
. 5^ act diagonalization (ED) level [8|; the LMG model can 
be solved by both MFT and ED. BECs offer the excit- 
ing possibility of realizing macroscopic superposition (MS) 
states (NOON states) with tens to thousands of parti- 
cles, and thereby pushing the limits of quantum mechan- 
ics [9]. However, the physical context of a BEC in a double 
well has quite different regimes. Thus different approaches 
are required, including MFT [TUHT2] . the multiconfigura- 
tional time-dependent Hartree (MCTDH) theory [LUfTj] . 
and our own [T5] and others' [TB1[T7] work on ED of the 
LMG model with additional terms, including tilt or bias. 
Without the active use of tilt, MS states are destroyed by 
potential fluctuations [TU1IT5] . 



In this Letter, we use a biased LMG model to investigate 
the quantum sloshing of many bosons in a tilted double- 
well, with the goal of guiding dynamical creation of MS 
states in BEC double- well experiments. Quantum sloshing 
is the tunneling dynamics of a system in which all atoms 
are initially localized in one well. At the quarter and three- 
quarter periods of the ensuing cyclical dynamics, one finds 
an MS state. The basic concept of quantum sloshing is 
similar to what occurs in an rf SQUID, in which the many 
body wavefunction oscillates between two macroscopically 
distinct states [18]. 

The LMG model is applicable for 



X=l(N 2 ~l)U]/(2hiu)<l 



(1) 



1 Sometimes called "two-mode Bose Hubbard" or "two-site Bose 
Hubbard" in cold quantum gases literature, although the LMG 
model precedes the Bosc-Hubbard model by 24 years. 



where A is the number of atoms, U is the interaction en- 
ergy, and u> is the local trap frequency in each of the two 
wells [inj[TH]. Criterion ([1]) means that only the lowest 
single-particle state in either well is occupied. Addition- 
ally, it is required that all dynamical perturbation have an 
energy much less than fvjj. Under these restrictions, there 
are two regimes: £/A = J/N\U\ 3> 1, where J is the tun- 
neling energy, we denote the Josephson regime; ( < 1 we 
denote the Fock regime [5]. It is in the latter that MS 
states occur. The experimentally observed phenomenon 
of self-trapping, in which the system becomes stuck in 
one well, can be attributed to long tunneling times [20] . 
These oscillations with long tunneling times are associated 
with the presence of Schrodinger-cat-like, or NOON-like, 
MS states. Such states collapse to a narrow distribution 
of Fock states in the presence of a very small tilt fT5] . 
and tunneling is exponentially suppressed. Past studies 
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on the suppression of tunneling have focused on environ- 
mental effects such as finite temperature or coupling to 
a reservoir [2"Trf2"3] . Therefore, in the presence of a small 
tilt, finite temperature, or coupling to a reservoir, the MS 
states are dynamically unaccessible. Tilt displays radi- 
cally different behavior than these other forms of suppres- 
sion of the tunneling. Namely, tunneling resonances occur 
for critical values of the tilt when the barrier is high. At 
resonance, MS states states reappear [15], and therefore 
tunneling is again observed. We show that at tunneling 
resonances the oscillation time between wells is hundreds 
of orders of magnitude faster and less sensitive to devia- 
tions in the tilt than in the symmetric case. 

This speed-up permits us to propose a simple scheme 
for the creation of MS states both for few-body and many- 
body systems. Whereas past proposals involved ramping 
the barrier height [24] or continuous variation of atom- 
atom interactions via Feshbach resonance [22], MS states 
are realized periodically in our scheme when all param- 
eters are fixed. These MS states take the form of pro- 
tected NOON states. Using the tunneling times of quan- 
tum sloshing to positively identify NOON or NOON-like 
states has been a significant tool in the study of SQUID- 
based NOON states, where the double well consists of 
right- and left-circulating states on a ring [18] . In cold 
atom experiments, quantum tunneling effects have so far 
been restricted to one or two atoms [2"6ll2"7] , and MS states 
of NOON-like form have yet to be clearly identified. 

We briefly mention methods for regimes other than 
that of the biased LMG. For \ > 1, MFT is applica- 
ble for C/N ^ 1; in this regime, called the "linear" or 
"Josephson" regime the double-well system can be 
described by a pendulum in a 2D phase-number phase 
space [MI21I21I2E1I2H]- In contrast, for (/N > 1, called 
the "nonlinear" or "self-trapping" regime, MFT methods 
find macroscopic self trapping in one well [T0UT21I29] . The 
concept of phase is well-defined in the Josephson regime, 
together with a clear scmiclassical limit; it is not well- 
defined in the "nonlinear" regime [121 ED], and then re- 
quires a more strongly quantum approach |31U32j . 

MCTDH theory is such an approach. It does indeed 
find different dynamics from mean-field predictions, par- 
ticularly in the Fock regime |33) . This purely numeri- 
cal method is more exact [2] than MFT or LMG ap- 
proaches [U, 20, 24 ( 25, 34,(35] and in principle superior, 
as it can span all regimes. However, our modified LMG 
approach has the advantage that it is simpler and leads 
to exact and perturbative analytical expressions together 
with straightforward simulations; a mix of analytical and 
numerical methods at different levels of approximation is 
useful. Moreover, because our methods are not computa- 
tionally intensive they can be easily extended to two and 
three dimensions |19| , unlike the more exact but computa- 
tionally demanding MCTDH theory Under condition ([T]), 



2 Sols and Lcggett subdivide this regime into "Rabi" and "Joseph- 
son" ; we use only the term "Josephson" for simplicity. 



MCTDH and LMG methods should converge. 

Thus we work with the biased LMG model. Although 
mean-field theory for the asymmetric trap has been con- 
sidered before [11], the MFT approximation is, as we have 
said, a poor one in the Fock regime, £ < 1. Experimen- 
tally, tilt appears both as a systematic error and deliber- 
ately in device applications [2"6ll2"Tll3l)rl38j . Tilted optical 
lattices are especially relevant to applications in gravit- 
ometry [36], quantum computing |26il39j . and atomtron- 
ics [55] . 

The two-mode Hamiltonian for N weakly interacting 
bosons in a tilted two-well potential, or biased LMG 
model, is 

k = - J E %' +Uj2™i (»i - !) + AVh ^ ( 2 ) 

JW 3 

where the subscript j € {L, R} is the well or site index, 
J is the hopping strength, U is the interaction poten- 
tial, and AV is the tilt. Here bj and satisfy the usual 
bosonic annihilation and creation commutation relations 
and hj = bjbj. Eq. ([2]) can be derived from first-principles 
quantum field theory for weakly interacting bosons at zero 
temperature. An arbitrary state vector in Fock space is 
given by 

= En i= °n L \ n L, N- n L ) , (3) 

where is the number of particles in the left well and 
c 1lL £ C. We require the total number of particles N to 
be constant. Under this restriction, the Hamiltonian re- 
duces to an (N + 1) x (N + 1) tridiagonal matrix [24] . 
We consider the dynamics of a system in which all par- 
ticles initially occupy the right well, i.e., \ip(t = 0)) = 
|0, N). In the Schrodingcr picture, the time evolved kct 
is \ip(t)) = exp(-iHt/h)\ip). The probability of find- 
ing riL particles in the left well at some time t > is 
P nL (t) = \{n>L, N — n^|?/'(t))| 2 , the average occupation of 
the left well is nj,(t) = (ip(t)\hL\'ip(t)}, and the average 
variance is o- 2 L (t) = {ip{t)\h 2 L \ip{t)) — n 2 L . 

We first consider the simple case of noninteracting par- 
ticles, U = 0, in a symmetric potential, AV = 0, to il- 
lustrate the problem. This case is exactly solvable. The 
probability of finding all particles in the right well, i.e., 
rii = 0, is 

P (t) = cos 2W ( Jt/ft). (4) 

The tunneling period is T = wh/J, which is independent 
of JV. When t = T/2, the system is in state |7V, 0) and 
all particles have tunneled into the left well. The average 
occupation and variance of the left well are 

n L (t) = Nsm 2 (Jt/h), (5) 
= (N/4) S m 2 (2Jt/h). (6) 

The particles therefore tunnel sinusoidally between wells 
with a frequency 2J/h. The variance is greatest when t = 
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Fig. 1: (Color online) Suppression of tunneling for noninter- 
acting atoms, (a) Shown are the probability densities P„ L (t) 
(colorbar) for all number states when N = 100, U = 0, and 
AV — 2.1 . Only N/2 = 50 particles tunnel between wells, (b) 
The tunneling amplitude and (c) the frequency of oscillation as 
a function of tilt. When AV > 2JV N — 1, tunneling is com- 
pletely suppressed. Particles tunnel between wells faster in a 
tilted potential than in a symmetric potential. 



T/4. At this time, the probability of finding nr, particles 
in the left well is 



Pnr.iT/ A) 



-N 



N\/[n L \(N - n L )\] 



(7) 



The system is in a truncated coherent state, i.e., a binomial 
superposition of all number-states. 

The tilted case is also straightforward [H], although 
to our knowledge the following expressions have not yet 
appeared in the literature. The dynamics are slightly dif- 
ferent when AV ^ 0, since the occupation of the left well 
now is 

n L (t) = Asin 2 (ujt/2), (8) 



where the amplitude and frequency of oscillation are 
A = N/[1 + (AV/2J)% 



lo = (2J/H) v / l + (AV/2jy 



(9) 
(10) 



When = 2 J, only N/2 particles tunnel between wells. 
Fig. QJa) shows the probability densities P nL (t) in this 
case. In Figs. QTb) andCTTc), Eqs. © and §TU§ are plot- 
ted as a function of AV. From the expression for A we 
observe that tunneling between wells is completely sup- 
pressed when |AV| > 2J^/N — 1. Because the hopping 
strength J is much smaller than the barrier height, tun- 
neling is highly sensitive to small tilt. 

We proceed to consider how a small interaction term 
changes this scenario in the Josephson regime, £ 3> 1, in 
a symmetric potential, AV = 0. For the non-interacting 
system a single frequency 2 J/h characterizes n l (t) ; in con- 
trast, N dominant frequencies emerge in the interacting 



Fig. 2: (Color online) Damped tunneling in the Josephson 
regime, (a) Probability densities P nL (t) for all number states 
when TV = 10 and (/N = J/NU = 10 for t < T 1/2 . All 
particles tunnel between wells with period T — nh/J. (b),(c) 
Average occupation (top panel) and number variance (bottom 
panel) of the left well for longer times, (b) Oscillations between 
wells are damped by atom-atom interactions, (c) The first tun- 
neling revival occurs when t = T r = nh/U. The colorbar is the 
same as in Fig. [TJ 



system in the Josephson regime. The average occupation 
of the left well is given by the modulated signal 

n L (t) = (N/2) [1 - cos(2Jt/h) cos^iUt/K)] , (11) 

to lowest order in perturbation theory [T§] in N/£. We 
have also verified this result through simulations. In 
Eq. (fTTj) the high frequency carrier depends only on the 
hopping strength J while the low frequency envelope de- 
pends on both the interaction potential U and the total 
number of particles N. The envelope reaches half its max- 
imum value when 



t = T 1/2 = (ft/(7)cos- 1 [2- 1/(Ar - 1) ] 



(12) 



At times t <C I\/ 2 , all particles tunnel between wells with 
period T, as in Fig. HJa). At times near T1/2, on the 
other hand, only half the particles tunnel between wells 
with period T. When t ~ 2T , 1 / 2 , there is essentially no 
tunneling (see Fig. [2Kb)). Small interactions thus damp 
the oscillations between wells [201140]. However, tunneling 
revivals occur periodically with period T r = Tth/U. The 
first tunneling revival occurs when |i — T r \ < Ti/ 2 , as 
shown in Fig. ^c). The separation of time scales, T X j 2 <C 
T r , occurs only for N » 1, as evident in Eq. (fT2"j) . 

For the remainder of our discussion, we turn to the high 
barrier limit, £ <C 1, as it is key to the dynamic production 
of MS states. We assume U > without loss of generality 
with respect to the dynamics. Using perturbation theory, 
it can be shown [T51[P9~| that the eigenstates are MS states 
of the form \<j>±) = (\N - n L , n L ) ± \n L , N - n L )) /y/2 
to lowest order in £. The degenerate number states in 
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the f = limit split into symmetric and antisymmet- 
ric MS states for small £, with an energy difference of 
AEjy- 7lL . The two eigenstates with the highest eigen- 
value are nearly-degenerate MS states of the form \<fi±) = 
(| AT, 0) ± |0, N)) /a/2. The energy difference between \<j)±) 
is 

AE N = 4[/(C/2) JV A/[(A- 1)!]. (13) 

The characteristic frequency is wjv = AEn/H. Notice that 
since AEn is a very small number (£ <C 1), lon is also very 
small, and decreases rapidly with increasing N. 
All particles occupy the right well with probability 



P (*) = 1 - Pw(t) = cos 2 {oj N t/2), 



(14) 



In Fig. EJa), we plot the probability densities P nL {t) and 
the average occupation nr,(t) as a function of time. The 
tunneling period is Tn = 2-k/lun. The average occupation 
and variance are 



n L (t) = N sin 2 (uj N t/ 2), 
r 2 n Jt) = (N 2 /4)sm 2 (u; N t). 



(15) 
(16) 



In this regime, as in the noninteracting case, all N par- 
ticles oscillate sinusoidally between wells. There are two 
important differences. The first is that the period of os- 
cillation depends on N and can become quite large for 
large values of N. Note that, for large N, self-trapping is 
observed for exponentially long times, in agreement with 
mean-field approaches. Second, at time t = T/v/4, we find 
that P/v = P = 1/2. At this time, all particles simulta- 
neously occupy both wells and the system is described by 
a NOON-like MS state, a new prediction that cannot be 
achieved in the mean-field limit, and which is intriguing 
in the framework of quantum information. 

Guided by this interest, let us characterize the entan- 
glement at t = T/v/4, by utilizing four standard entan- 
glement measures: the average local impurity, or "Q- 
measure" [H] , the local entropy of entanglement [55] , the 
Schmidt rank k |43| . and a "macroscopic superposition 
size" (MSS) based on physical measurement [H] . The Q- 
mcasurc is given by Q = [(N+l)/N][l-(Trp 2 L +Tip 2 R )/2], 
where p L ( R) = Ti R ^\ip(t))(tp(t)\. The entropy is S = 
~Yl nij =o P n L (t)^og N+1 P nL (t), and tnc Schmidt rank k 
is given by the number of non-zero eigenvalues of the re- 
duced density matrix pl- Finally, the MSS measure is 
Cs = N/n m i n , where n m in is the minimum number of par- 
ticles one must to measure to distinguish both branches, 
i.e., \N — til, ni) and [ul, N — Ul}- These measures take 
the value Q = S = 0, k = 1, and Cs — if and only if 
\tp(t)) is a pure state. This occurs at t = Tv/2 and Tv. 
At time T/v/4, we find that each measure reaches a max- 
imum value of Q = N/[2(N + 1)], S = log N+1 (2), k = 2, 
and C 5 - N/{n L + 1). 

We proceed to consider the effects of tilt or bias in the 
LMG model. For ( <g; 1 the tunneling between wells is ex- 
tremely sensitive to tilt AV. Using perturbation theory, 
it was shown that when AV > 2AEn/N the eigenstates 
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Fig. 3: (Color online) Tunneling resonances in a few-atom sys- 
tem. Shown are the probability densities P„ L (t) when N — 7 
and C = 0.1, for (a) AV = and (b) AV = 4(7. (a) All particles 
tunnel between wells with period Tn- At time t = Tat/4, the 
system is described by an MS state, (b) Only N — 2 — 5 parti- 
cles tunnel between wells. The oscillation frequency is 5 orders 
of magnitude faster than the symmetric case, (c) Tunneling 
amplitude as a function of tilt AV for N = 5 and ( = 0.1. 
Tunneling resonances occur when AV = AV P = 2pU. At res- 
onance, N — p particles tunnel between wells. The insert is a 
zoom around AV/2U = 2. 



are nearly perfect number-states of the form \nr,, N — nr) 
and \N—riL, nr). In this case, since the number states are 
near-eigenstates of the Hamiltonian, the initial condition 
\tp{0)} is nearly stationary and tunneling between wells is 
strongly suppressed. This is quite different from suppres- 
sion of tunneling due to thermal effects or coupling to a 
reservoir [21H23| : our system is closed and suppression is 
due to an internal parameter, namely, imperfections in the 
trapping potential. 

However, for stronger tilts quasi-degenerate MS eigen- 
states of the form |</>±;p) = (|JV - p, p) ± |0, N)) /y/2 
reappear [T3| for AV = AV P , with a splitting in the ener- 
gies between the states \4>±;p) ec 



equal to AE P N: where 



AE% 



4U((/2) n ~p{N -p) 
{N-p-l)\ 



N\ 



p\{N -p)\ 



(17) 



In this case, the potential difference is exactly compen- 
sated by the repulsive interaction of p particles in the 
lower well. Then, at resonance, the tunneling frequency is 
AE^/h. The average occupation of the left well is 



N 



n L (t) = (N -p)sm 2 {uj p N t/2) 



(18) 



to lowest order in £. Here, N — p particles tunnel between 
wells with period = 2tt/w^. At time t = T%/2, N - p 
particles are in the left well. To compensate the tilt, p 
particles remain in the right well at all times. When t = 
Tjy/4, the system is described by an MS state such that 
Pp = Pa = 1/2 and the entanglement measures Q, S, and 
k reach the same values as in the symmetric case. 
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Fig. 4: (Color online) Tunneling periods in a many-body sys- 
tem, (a) Tunneling period Tjv versus the total number of par- 
ticles N when £ = 0.1 and AV — 0. For large N, tunneling 
becomes very slow, (b) At resonance, AV = AV P = 2pU, only 
N — p particles tunnel between wells. Shown are the tunneling 
periods versus p for N = 40 to 100 with £ = 0.1. At reso- 
nance, the oscillations can be hundreds of orders of magnitude 
faster than in the symmetric case. 



The entanglement measured by Cs is again C$ = N, 
provided that the measure is taken in the higher well. In- 
terestingly enough, if the measure is taken in the lower 
one, the entanglement measured is smaller, Cs = N/(p + 
1), since p + 1 particles must be measured to distinguish 
between both branches. In Fig. [3Jb) the tunneling dynam- 
ics for the second resonance, i.e., p = 2, are illustrated for 
a system of N = 7 particles. Near a resonance, tunneling 
is suppressed when 



\AV-AV p \>2AE p N /{N-p). 



(19) 



This is due to the fact that, near the resonance, the eigen- 
states are again near-perfect number states when the devi- 
ations of the tilt with respect to the resonant one exceed 
this quantity [TS]. Then, the condition for the suppres- 
sion of the tunneling in the symmetric case is that the 
tilt exceeds 2AE^/N while in the case of the asymmetric 
potential, when the tilt coincides with a resonance, the 
condition for the suppression of the tunneling is that the 
difference between the tilt and that of the resonance ex- 
ceeds 2AE P N /(N -p). 

Moreover, tunneling near resonance is much faster than 
tunneling in a symmetric potential since both tunneling 
frequencies, uj p n and uin, are proportional to AE p N /h and 
AEn/H respectively, and we have shown that AE 1 ^ is 
greater than AEn by many orders of magnitude. In 
Fig. |4ja), we show the symmetric tunneling period Tjv 
versus N when £ = 0.1. Clearly, TV becomes very long 
as N becomes large. For instance, in a typical symmet- 
ric double- well used in experiments [55], 200 87 Rb atoms 
tunnel between wells with period T200 = 1-15 x 10 635 ms 
when ( = 0.0964. Furthermore, tunneling is completely 



suppressed for deviations in the tilt greater than 4.16 x 
10 -636 nK • fee. Obviously, one does not expect to observe 
many-body tunneling in this regime. Notice that, under 
the same conditions, systems with as few as N = 1, 2, and 
3 87 Rb atoms yield tunneling times as long as T\ = 466 ms, 
T 2 = 4840 ms, and T 3 = 134000 ms, respectively. Even in 
a few-particle system, tunneling times can be prohibitively 
long, demonstrating self-trapping behavior @]. 

However, tunneling at resonance can be hundreds of or- 
ders of magnitude faster than the symmetric case, as in 
Fig. |U[b). For the 200-atom system discussed above, when 
p = 197, we find that N — p = 3 particles tunnel between 
wells with period T^qq = 117 ms. This resonance occurs 
when AV = AV lg7 = 210 nK • k B - An MS state of the 
form = (|3, 197) - i|0, 200))/V2 will be realized at 
r 20o7 4 = 29 - 25 ms - Likewise, we find that Tj 98 = 34.3 ms 
and T^ 9 = 33.0 ms when AV = AV W& = 211 nK ■ k B and 
AV499 = 212 nK • k B , respectively. Finally, at resonance, 
this system is sensitive to deviations in the tilt on the order 
of 0.273 nK-k B , 1-40 nK-fc B , and 2.90 nK-fc B forp= 197, 
198, and 199, respectively. Thus, the observation of the 
tunneling of a few Rb atoms is made possible by tunnel- 
ing resonances in a many-body system. Moreover, embed- 
ding NOON-likc MS states in the many body wavefunction 
leads to the possibility of larger MS states of more than 
three particles; the scaling of the advantage in tunneling 
time gained via tilt can be calculated as r = AE^ / AE P N , 
from Eqs. (fl~3"|) and (fT?]) . Taking p' = N — p as the number 
of atoms in the embedded MS state, and using Stirling's 
approximation Inn! ~ nlnn — n. we find 



In: 



(- Inn + 1 + In §)n + + ^ - \ In n)p' 



+(|lnp'-| + iln4-lnCy 



(20) 



the second set of cross terms is vital. This expansion 
matches the exact expression very well, as illustrated in 
Fig El 

In conclusion, we used the two-mode approximation to 
develop a Fock space picture of a system of ultracold 
bosons in a tilted two-well potential, covering all regimes 
of barrier height, from the Joscphson regime to the regime 
in which one can find NOON states. In the latter regime, 
which occurs when the barrier is high, a small tilt causes 
the complete suppression of tunneling, leading to self- 
trapping. Long tunneling times prevent the observation of 
many-body tunneling even in a symmetric potential and 
MS states are too sensitive to fluctuations in the trapping- 
potential to be realistic. However, in this regime, tunnel- 
ing resonances occur when the tilt can be compensated by 
atom-atom interactions. At resonance, tunneling is much 
faster and less sensitive to tilt than in a symmetric po- 
tential. Furthermore, tunneling resonances can be used to 
create NOON-likc MS states embedded in and protected 
by a larger many-body system. 
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Fig. 5: Embedding the MS state in the many body wavefunction. 
Relative increase r in tunneling time and robustness against 
potential fluctuations, for an MS state of N — p — 10, 50, 100 
particles (left to right); plotted as a function of the number of 
particles in the full state, N, for exact (solid black curves) and 
approximate (dashed red curves) expressions, all for £ = 0.1. 
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